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Abstract 

The equations of state of spin-polarized nuclear matter and pure neutron matter are 
studied in the framework of the Brueckner-Hartree-Fock theory including a three-body force. 
The energy per nucleon Ea{§) calculated in the full range of spin polarization 5 = P] p Pl for 
symmetric nuclear matter and pure neutron matter fulfills a parabolic law. In both cases 
the spin-symmetry energy is calculated as a function of the baryonic density along with the 
related quantities such as the magnetic susceptibility and the Landau parameter Gq. The 
main effect of the three-body force is to strongly reduce the degenerate Fermi gas magnetic 
susceptibility even more than the value with only two body force. The EOS is monotonically 
increasing with the density for all spin-aligned configurations studied here so that no any 
signature is found for a spontaneous transition to a ferromagnetic state. 



1 



I. INTRODUCTION 



Studies of spin-polarized nuclear and neutron matter have been mainly focussed on the 
possible onset of a ferromagnetic transition in the neutron star core. This transition could 
explain in fact the high intensity magnetic fields (10 12 gauss) estimated from the timing 
observations in pulsars and magnetars (for a review see Ref. [||]). Besides this exciting issue 
the motivation for such studies can be based on a more general context of nuclear physics. 

First of all, from the change of the energy per nucleon with the spin polarization one may 
extract a theoretical prediction for the spin symmetry energy, whose empirical value is so 
far quite uncertain. The nuclear matter instability against spin fluctuations is driven by the 
Landau parameter Go which is determined from the spin-symmetry energy. The value of this 
parameter is still a largely controversial topic and no agreement exists among the different 
approaches to the theory of nuclear matter |2j. Experimental information could come from 
spin giant resonances, which have not yet been clearly observed. Information could also 
come from heavy-ion collisions as soon as polarized heavy targets become available. 

The second issue related to the study of spin-aligned states of nuclear matter is the prop- 
agation of neutrinos in neutron stars. It has been shown that the neutrino mean free path 
is strongly affected by the magnetic susceptibility. The latter is sizeably suppressed by the 
strong correlations in nuclear matter and, as a consequence, the mean free path might change 
sizeably and, eventually drop to zero in the presence of a ferromagnetic transition ||||] . 

There is a guess that the ferromagnetic transition could be a relativistic effect due to 
7r-exchange and in fact all calculations, based on the relativistic mean field approch, pre- 
dict this transition to occur in dense matter f||(|. On the other hand, non- relativistic 
approaches |]7HlO| do not support such a transition except Hartree-Fock calculations with 
phenomenological Skyrme-like forces (for a review see Ref. ||). This aspect cannot be dis- 
connected from the problem of the in-medium nucleon-nucleon (NN) force, which is poorly 
known in dense matter due to the lacking of empirical constraints far above the satura- 
tion density. However important relativistic effects can be incorporated into the effective 
interaction via the three-body force associated with a virtual nucleon-antinucleon excita- 
tion jnj. Moreover non- relativistic calculations including only two-body forces miss the 
empirical saturation point of nuclear matter [12|. So it seems worthwhile to investigate the 
spin-aligned states of nuclear and neutron matter in the non-relativistic Brueckner theory 
with three-body forces. They contain not only the above mentioned relativistic contributions 
but also nucleonic excitations which decisively enhance the agreement between theoretical 
and empirical saturation density |13| . 
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II. FORMALISM 



The spin and isospin asymmetric nuclear matter (ANM) consists of spin-up neutrons 
(n |), spin-down neutrons (n |), spin- up protons (p f) and spin-down protons (p I) in 
different density states: p n |, p n ^ p p ^, and p p ^ respectively. Therefore four parameters are 
required to specify a given configuration of spin and isospin ANM. The Fermi momenta of 
the four components are generally different from each other, and related to their respective 
densities by the following relation: 

where A denotes the z-components of isospin and spin, i.e., A = (t z , a z ). Instead of p\, one 
can use the following four parameters to identify a given spin and isospin state, 

n _ Pn - Pp r _ Pn\ ~ Pnj r _ Py\ ~ Ppj 
P ; On ; Op , 

P Pn Pp 

where p, p n , and p p are total density, neutron density and proton density, respectively. The 
ratio (3 is the isospin asymmetry parameter and 5 n and 5 P are the spin asymmetry parameters 
for neutrons and protons, respectively. 

The starting point of the Brueckner-Bethe-Goldstone (BBG) approach is the reaction 
G-matrix. The G-matrix incoporates strong short-range correlations in nuclear medium by 
means of the infinite ladder diagram summation of the bare NN interaction. It satisfies the 
Bethe-Goldstone equation. The latter can be expressed for the spin-isospin ANM in the 
total angular-momentum basis as follows, 

G T L S J> XX \u, P; q, q'; p, (3, 5 n , 5 P ) = v T L S J{q, + / g" 2 dg"^(g, q") 

H T „ J 



(Q XX (q",P)} n TSJ,X,\' ( r> _// Joe X \ 

X «,-e$(q»,P)+ir, Gl "» ^^^^ P ' & ^ 

where vj i i/(q, q') are the partial wave components of the NN interaction, P = ki+k 2 = k'i+k' 2 

ki — k 2 k' — k! 

is the total momentum, q = and q ' = —^—^ — - the relative momenta of the two 

particles in their initial state and final state, respectively. The Pauli operator Q xx '(q",P) 
and the energy denominator e x 2 (q", P) = eX (k'{) +t x ' (k' 2 r ) have been angle averaged in order 
to remove the coupling between different channels a = {J ST}. It is worth noticing that the 
different components of the G-matrix differ in general from each other due to the dependence 
of the Pauli operator and energy denominator on the spin-isospin configuration (A, A '). The 
single particle energy is given by e x (k) = h 2 k 2 /2m + U x (k). The continuous choice for the 
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auxilary potential U x (k) is adopted in the present calculations since, on the one hand, it has 
been shown to provide a much faster convergence of the hole-line expansion than the gap 



choice [TE], on the other hand, it decribes physically the single-particle potential felt by a 
nucleon in nuclear medium. In the continuous choice, U x (k) is the real part of the on-shell 
mass operator, i.e., 

U\k) =ReJ2 n x '{k')(k\, k 'X'\G[u = e x {k) + e x '(k'), P]\kX, k'X') A , 
k',\> 

where the subscript A denotes antisymmetrization. For spin and isospin ANM, it is conve- 
nient to split it into two contributions as 

U x = U xx + U xx ' , (A'/A). 

Each individul contribution is calculated by casting it into the partial wave expansion, 
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q,q;p,P,8 n ,8 p ). 

The summation over partial wave states is physically constrained by the selection rule S + 
T + L = odd due to the Pauli principle and consequently the antisymmetrization of the G- 
matrix simply implies multiplication by a factor of 2 for the allowed partial wave channels. 
For spin symmetric case (8 n = 8 P = 0), a spin-up neutron (proton) has the same Fermi 
momentum as a spin-down neutron (proton) and thus the single-particle potential felt by 
a nucleon does not depend on the direction of its spin. The summation on the spins of 
the two particles in the final state and the average of that in the initial state remove the 
non-diagonal contributions in angular-momentum from the single-particle potential. One 
easily finds 

U^<(k) = - E U^'^ik) 

a z a' z 
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For spin-asymmetric but isospin-symmetric nuclear matter, we have (3 = and 8 n = 5 P 
Pl~ Pi 



P 



8. In this case, the single particle potential becomes. 
2T + 1 



EE 



■L-L' 



C{—a z —a'\SS z 
y 2 2 1 



2Gl S L ?' azCr '*(uj,q,q;p,8) 



TSJ " LL' S z 

x Y, C ( L ' M lSS z \JM l + S Z )C(LM L SS Z \JM L + S z )Y* ML (q)Y LML (q). 



4 



The present calculations will mainly consider this spin-polarized nuclear matter as well as 
the spin polarized neutron matter. 



III. RESULTS AND CONCLUSIONS 

We performed some calculations within the BHF self-consistent approach above de- 
scribed. The Argonne Vis force is adopted as bare two-body interaction. This has been 



implemented by a microscopic three-body force , which is described in detail in Ref. [16] to- 
gether with the average procedure to transform it into an effective two body force. In Fig. 1 
the energy shift per nucleon E A {S,p) — E A {0,p) in symmetric nuclear matter is reported 
as a function of the square of spin polarization S 2 for a set of densities. Due to the linear 
dependence on <5 2 , also reported in Refs. PJIOfl one can write: 

E A {5,p) =E A {0,p) +S sym {p)5\ 

i.e. the spin dependence of the energy per nucleon can be simply expressed in terms of 
spin-symmetry energy 

l d 2 E A (8,p) 

c sym\P) 2 Qp 

in the density range here considered. 

The <5 2 -law is mainly due to the BHF approximation (two-hole line terms only). The same 
behavior is in fact exhibited by the energy vs. isospin either []14||. The effects of three- hole 



line terms are rather small when adopting the continuous choice for the auxiliary potential 



lq| ; this choice is also adopted for the present calculations. The slope of the energy shift 
is monotonically increasing with density so that no signature for a ferromagnetic phase 
transition in symmetric nuclear matter is expected. The effect of three-body force is to 
enhance this slope for densities above the saturation point. This effect is more clearly shown 
in the plot of £ sym vs. density in Fig. 1. 

The magnetic susceptibility has been also calculated from the spin-symmetry energy 

p 2 p 

X 



2£ 

where p, is the average of neutron and proton magnetic moments (in neutron matter p, is 
the exact neutron magnetic moment). Usually one calculates the ratio of x to xf, Xf being 
the magnetic susceptibility for a degenerate free Fermi gas. The effect of strong correlations 
in nuclear matter due to the two-body force is a reduction of x with respect to xf- This 
reduction increases with density up to a factor of 0.3 at p = 0.8/m -3 . The above result is 
common to most Bruckner calculations [§,0J||1O|. More pronounced is the queching due to 
three-body force. 
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Fig. 2 also shows the Landau parameter Go describing the spin density fluctuations in 
the effective interaction. Go is simply related to the spin-symmetry energy or, equally, to 
the magnetic susceptibility by the relation 

X ™* 



Xf 1 + G 

where m* is the effective mass. A magnetic instability would require Go < — 1 which is 
analogous to the condition F < — 1 for the mechanical instability giving rise to the liquid- 
vapour phase transition. But, the value of Go vs. density from the BHF calculations is 
always positive and monotonically increasing up to the highest density. The three-body 
force pushes up the curve of \. This result is in strong disagreement with the prediction 
with Skyrme forces. This is not a complete surprise since Skyrme forces are only well suited 
in the proximity of the empirical saturation point. Astonishing is the strong disagreement on 
this respect with the relativistic approaches because the three-body forces contain already 



important relativistic effects ||16|| . The accurate knowledge of Go should lead to reliable 
predictions on the spin and spin-isospin giant modes as well as spin-spin part of the optical 
potential [|J. 

The above calculations have been also repeated for the case of pure neutron and reported 
in Fig. 2. The same conclusions can be drawn as to the absence of the ferromagnetic phase 
transition and the quenching of the magnetic susceptibility caused by the strong correlations 
from the two- and three-body forces. This quenching should have a strong influence on the 
neutrino propagation in dense matter such as supernovae and neutron stars. In the case 
of the transition to a ferromagnetic state it has been shown that the mean free path could 
drop to zero || that could have remarkable consequences as, for instance, on the neutron 
star cooling. 
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FIGURES 



FIG. 1. Uper windows: The EOS of spin-asymmetric nuclear matter as a function of 
spin-asymmetry at five values of density, predicted by Brueckner Hartree-Fock calculations adopt- 
ing pure AVi8 two-body force (right-uper window) and AVis plus the TBF(left-uper window). 
Right-Lower window: density dependence of spin-symmetry energy for both cases with the TBF 
(solid curve) and without the TBF(dash curve). Left-lower window: Magnetic suscepbility x/xo 
(curves with symbol) and Landau parameter Go (curves without symbol) as functions of density. 



FIG. 2. The same as Fig.l for spin-polarized neutron matter. 
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